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ABSTRACT 

Extrapolations of numerical data obtained from finite cluster calcula- 
tions to the infinite volume limit can often give incorrect results. We discuss 
four separate cases: (a) the intensity of the lowest two-photon absorption in 
the infinite polyene, (b) bond alternation in the infinite polyene, (c) Cooper 
type pairing in the simple Hubbard model, and (d) pairing within the extended 
Hubbard model. 
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I. Introduction 

Approximate analytic techniques have been known to lead to qualita- 
tively incorrect predictions for solid state systems with strong electron corre- 
lation. Exact numerical techniques are therefore often employed to determine 
the behavior of such systems. This approach consists of two steps, the first 
of which involves an exact calculation for a finite cluster. The second step in- 
volves an extrapolation of the finite cluster data to the thermodynamic limit. 
We point out that extrapolations are not always valid, and in certain cases 
can lead to erroneous theoretical predictions. 

We believe that there is a strong necessity for such a discussion cur- 
rently. The technique of using exact numerical data for finite systems to 
determine the behavior in the solid began with the classic calculations for the 
Heisenberg chain by Bonner and Fisher 1 , who, however, emphasized the ne- 
cessity of having one dimension and short range interaction. In recent years, 
however, the numerical approach has been employed more and more liberally, 
and calculations are routinely being done for higher dimensions, as well as for 
Hamiltonians with long range interactions. Our hope is that the present work 
will lead to more systematic studies of finite size effects. 

We discuss four specific cases here. The correct solutions to these 
problems are of interest to chemists and physicists. These four different cases 
demonstrate different finite size effects, and in each case we show that literal 
interpretation of finite cluster calculations predict qualitatively incorrect result. 

II. Exact Finite Cluster Calculations 

A. The intensity of the two-photon absorption (TPA) to the 2A g and other 
low lying A g states in the infinite polyene The occurrence of the lowest two- 
photon state, the 2A g: below 2 the lowest one- photon state, the 1B U (note 
that optically inactive B u states need not be considered), in finite polyenes 
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necessarily requires a correlated description for the 7r-electrons 2 ' 3 . The entire 
discussion of TPA has, however, so far been limited to the energy of the 2A g 
state. Here we extend this discussion to the question of the intensity of the 
TPA to the 2A g (and other low lying A g states) as a function of the chain 
length N. 

This discussion will be in the context of the extended Hubbard Hamil- 

tonian, 

H = t £[1 - (-iyS](c+ a c i+h<7 + ct +1 ^ a ) 



1,(7 



+U £ n lA n hi + V ^(n, - l)(n i+1 - 1) (1) 

i i 

where all operators and parameters have their usual meanings. We begin with 
the U = V = limit. In Fig. 1 we show the occupancies of the highest two 
valence band levels and lowest two conduction band levels for the lA g (the 
ground state), the 1B U , the 2A g and the 2B U for arbitrary N. TPA is given 
by the third order nonlinearity x*- 3 ^— iv; u>, u, —ui) = Xtpai 

(3) {\A g \^\jB u ){jB u \ l i\kA g ){kA g \ii\lB u ){lB u \ii\lA g ) 
XTPA fa (" jBu ~ u) {u kAg ~ M (u lBu -u) U 

where [i is the dipole operator. Conventional wisdom has been that TPA 
occurs every time an energy denominator goes to zero in Eq.(2). We show 
below that this is not true. 

The intensity of the TPA to the 2A g was calculated exactly for U = 
V = and for N = 6, 8, 10 and 12, typical chain lengths for which calculations 
are done for nonzero U and V. As seen in Fig. 2., TPA intensity increases with 
N here, suggesting a large finite TPA intensity for — > oo. This, however, is 
incorrect. 

The model in Eq. (1) has charge-conjugation symmetry, which implies 
that for every single-particle valence band level at energy -e there exists a 



conduction band level at e. In the infinite chain, dipole allowed absorptions 
occur only between these symmetrically placed molecular orbitals, in agree- 
ment with ^-conservation. Thus for TPA to any intraband A g state, only two 
"symmetric" B u states need be considered. For the 2A g , these two B u states 
are the 1B U and the 2B U . For any N, the absolute values of the transition 
dipole moments (2A g \fi\lB u ) and (2A g \/j\2B u ) are equal, while because of the 
fermion character of the electron the products {lA g \/j,\lB u )(lB u \/j,\2A g ) and 
(lA g \n\2B u ) (2B u \n\2A g ) are of opposite signs. At N — > oo, these two prod- 
ucts become equal in magnitude, while having opposite signs. This implies 
that the contribution to in Eq.(2) by the 2A g vanishes 4 , and the TPA 
intensity to the 2A g and other band edge two-photon states is zero for N — > oo. 

We emphasize that the above analysis does not imply a zero intensity 
for the total TPA to intraband two-photon states. Even with the restriction 
to "symmetric" B u states, there can be very weak TPA due to relatively high 
energy A g states that are slightly removed from the band edge. For example, 
one can have such a high energy intraband A g state that is coupled to the 
1B U and a B u state that requires excitation from deep inside the valence band 
to deep inside the conduction band. Because of the small nonzero energy 
difference between the two B u states now, the cancellation is no longer total. 
Weak TPA to states away from the band edge is expected, and explains the 
shift of the TPA to higher energy from the optical band edge in the long 
chain 5 . Nevertheless, at infinite N, TPA to band edge two-photon states goes 
to zero, in contradiction to the prediction form Fig. 2. Since calculations of 
correlated chains are limited to the region where TPA increases with N, the 
integration of these two contradictory results is an important issue. 

For the noninteracting case TPA can be calculated for arbitrary N. 
In Fig. 3 we show the TPA intensity to 2A g as a function of N. It is seen 
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that the TPA intensity does go to zero at very large N, but only after a 
maximum is reached. Thus calculations near N = 10 predict qualitatively 
incorrect results. On the other hand, once the mechanism is understood, one 
can probe the N-dependence of where S~ and S + are the negative 

and positive contributions to the TPA. This is shown in Fig. 4, where the 
expected N-dependence is observed. 

We are now in a position to calculate TPA to 2A g in long correlated 
chains. We have adopted the novel configuration interaction (CI) approach 
of Srinivasan and Ramasesha 6 for the interacting Hamiltonian of Eq. (1), 
with a cutoff 4t in the energies of single-particle configurations. The approach 
can give the 2A g below the 1B U even in long chains. In Table 1 we compare 
calculated for many different U, V and 5. The ratio cannot be 1 in 

N = 20 - 30, but comparisons with the noninteracting cases strongly suggest 
that the TPA to the 2A g should be very weak in the long chain limit even 
for nonzero Coulomb interaction, in contrast to what would be predicted from 
exact evaluation 7 of the TPA near N = 10. This prediction is in agreement 
with the tiny TPA to 2A g in poly-BCMU polydiacetylene 8 , and the noted 
absence of the two- photon resonance in /3-carotene 9 . 
B. The bond-alternation problem in the one dimensional half-filled band 
Within the Su-Schrieffer-Heeger Hamiltonian, 

H = Y^[t-a(yi - y i+1 )}(c+ a c i+1 , a + cf +ha a, a ) + ^K^{y x - y l+1 ) 2 (3) 

i,u i 

bond-alternation in the infinite polyene is unconditional. Coulomb interac- 
tions can enhance this bond alternation. For long range Coloumb interactions 
H ee of the form , 

H ee = U ^2n^n hl + ^ j V j {n i - l)(n i+j - 1) (4) 

i i,3 
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two sharp inequalities give the ground state broken symmetry 10 . For off-site 
Coulomb parameters that obey the convexity condition, Vj+i+Vj-i > 2Vj, the 
ground state can only be a bond-order wave (BOW) with alternating bonds 
or a charge density wave (CDW) with a periodic modulation of the on-site 
charge density. The dominant broken symmetry is given by the inequality, 

\U + J2 V ^ (5) 

where a smaller left hand side (LHS) gives BOW, while a smaller right hand 
side (RHS) gives a CDW. All Pariser-Parr-Pople (PPP) parametrizations for 
the infinite polyene then predict enhanced BOW 10 . 

The approach involving extrapolation technique was also employed in 
this case, by Ramasesha and Soos 11 . These authors calculated the electronic 
energy gained for fixed bond alternation, when compared to the equal bond 
length case, for open even chains, closed N = 4n rings, and closed iV = An + 2 
rings, and extrapolated the energy gains in the three series against 1/N for 
the PPP-Ohno Hamiltonian. Enhanced energy gains compared to the simple 
Hiickel case was taken to be the proof of enhanced BOW. While the result 
for the carbon-based polymers agrees with the prediction from Eq. (5), let us 
examine this approach in detail. 

The extrapolation of the 4n-ring series is invalid, since the RHS of Eq. 
(5) in this case always has an extra term that precludes the CDW entirely and 
guarantees enhanced BOW. The even chain extrapolation has many different 
finite size effects. First, even chains have an odd number of bonds, thus 
precluding equal bond lengths for small N. Second, the range of the Coulomb 
interaction within PPP models is longer than the N for which calculations 
are being done. This leads to a different Hamiltonian for each N. Finally, 
real space many-electron configurations have entirely different matrix elements 
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of H ee within the open chain boundary condition and for the periodic ring. 
For example, the configuration ...20202020..., where the numbers denote site- 
occupancies, has a matrix element AU — 8V1 + 8V2 — 8V3 + AV4 in the ring, 
but the matrix element is AU - 1V X + 6V 2 - 5V 3 + AV 4 - 3V 5 + 2V 6 - V 7 in 
the open chain. This last difference shifts the boundary between the BOW 
and the CDW considerably, such that the open chain would predict enhanced 
BOW even as Eq. (5) predicts a CDW. Since the An + 2 ring series can only 
have a few data points, this makes the extrapolation technique itself suspect. 

We demonstrate the above by doing the open chain calculation for the 
PPP-Ohno parameters corresponding to polysilane. It has been assumed that 
this has a simple BOW ground state with very large alternation 12 . We choose 
a U of 9.04 eV and evaluate Vj from the PPP-Ohno relationship, 

U 

Vj = [1 + 0.3912 2 ] 1 / 2 (6) 

where the Rj are the distances between the polysilane "sites" , and are taken 
from literature 12 . Eq. (5) predicts a CDW, with different charge densities on 
the orbitals on the same silicon atom. The extrapolation technique neverthe- 
less predicts enhanced BOW, as seen in Fig. 5. 

Our motivation here was to demonstrate the inapplicability of the ex- 
trapolation technique for the bond alternation problem and not to claim that 
polysilane is a CDW. Whether or not polysilane is a CDW depends on whether 
the PPP-Ohno model is truly appropriate for polysilanes. It is, however, inter- 
esting to note that the occurrence of the lowest two-photon absorption consid- 
erably above the lowest one-photon absorption 13 occurs naturally within the 
CDW model. 

C. Pairing in Hubbard models A great deal of numerical effort has gone into 
investigating whether pairing of the Cooper type occurs in the weakly doped 
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two-dimensional Mott-Hubbard insulator. More recently, such a pairing mech- 
anism has been proposed 14 for doped Cqo. The two dimensional calculations 
have now been extended to the case of weakly coupled layers 15 . Specifically, 
these calculations involve determination of the ground state energies of the 
half-filled band system, and the electronic energies needed to add one and two 
electrons (holes). Pairing is supposed to occur if the quantity 

AE = 2E(N + 1) - E(N) - E(N + 2) (7) 

is positive. In spite of considerable effort, it is still not entirely clear whether 
pairing occurs in two dimension or not, or what the actual mechanism of 
pairing is in the various cases. 

We focus here on the mechanism of the pairing, and what this says 
about pairing in the infinite system. Fye et al 16 showed that pairing of the 
above kind is not limited to two or three dimension, but can also occur in one 
dimensional periodic rings with N = An atoms. The limitation to 4n rings is 
very suggestive, as small 4n-rings undergo Jahn- Teller distortions. In Fig. 6 
we show our results of the pairing calculations for N = 4 and 8. The pairing 
disappears for the Jahn- Teller distorted rings. To understand this better, we 
calculated the effect of the Hubbard U on the Jahn- Teller distortion of the 4n 
and 4n — 1 electron systems. These results are shown in Fig. 7. It is seen that 
U strongly reduces the gain in energy for the Jahn- Teller distortion of the even 
(4n) electron system, but the distortion of the odd numbered (4n -1) electron 
system remains virtually unaffected. This then gives us the mechanism of 
the pairing. The information regarding the tendency to distort is built into 
the wave function in the form of the various bond order parameters. By not 
letting the system distort, we are overestimating the energies in each case, but 
as Fig. 7 indicates, this overestimation is much stronger for the odd electron 
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case than for the even electron case. The calculated energy required to add 
one hole (or electron) to the undistorted system is therefore too high, thereby 
giving "pairing". 

This same mechanism can explain "pairing" in higher dimension. In 
Fig. 8 we show the results of our calculation for the simple cube, and for var- 
ious distorted forms. For the undistorted simple cube, we merely reproduce 
the results of White et al. 17 , — two distinct regions are found in which pairing 
is seen to occur. For the cube stretched along one axis, we see that the pairing 
still exists. This is not the Jahn- Teller mode for the cube, the noninteracting 
limit still having degeneracies. However, the pairing at small U has disap- 
peared for both the Jahn- Teller distorted cube as well as the structure with 
all three hopping integrals different. This clearly indicates the relationship be- 
tween pairing and the degeneracies in the single particle limit. We have once 
again calculated the effect of the correlation parameter U on the Jahn- Teller 
distortions. The result is identical to that in Fig. 7, the tendency to distort 
is almost unaffected for the 7 electron system, but for the 6 electron system 
this is strongly suppressed. The mechanism of pairing at small U in the cube 
is then the same as in the 4n ring. In the absence of distortion, the energy to 
add one hole or electron is being overestimated. 

We have repeated such calculations for a number of three dimensional 
molecules, and have found the following. For an N site system, pairing at 
small U occurs only if the (N — 1) (or (N + 1)) electron system has a strong 
tendency to have a Jahn- Teller distortion. Thus the mechanism of pairing at 
small U in all these case is related to the suppression of Jahn- Teller distortion. 
Since the Jahn- Teller distortion results from the discrete level degeneracies in 
finite molecules, we conclude that the observed pairing in in the small U region 
in similar calculations is a finite size effect. 
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Fig. 8 explains the oscillation in the pairing energy of the cube. Pairing 
at large U is unrelated to the tendency to have Jahn- Teller distortion. We have 
occasionally found such pairing at large U also in other systems, though not 
necessarily accompanied by pairing at small U. The occurrence of such large U 
pairing in nonbipartite systems (for e.g., in a pentagonal prism) indicates that 
this is also unrelated to antiferromagnetism. Currently we are investigating 
whether the large U pairing is also a finite size effect, and if so, the nature of 
the finite size effect. The present work indicates the absence of pairing in the 
two dimensional Hubbard model, as well as in weakly coupled layers. 
D. Pairing due to intersite interactions Many calculations exist that find a 
pairing driven by intersite interactions rather than the on-site Coulomb inter- 
action. Since the intersite interactions in these calculations are different, it is 
impossible to give a general critique. Nevertheless, we believe numerical cal- 
culations here would in most cases find phase segregation, rather than simple 
pairing, if adequate numbers of electrons or holes are added to the half-filled 
band. In the case of the t-J model, this tendency to phase segregate has been 
discussed 18 . We choose one simple example, a calculation claiming pairing 
driven by the nearest neighbor Coulomb interaction 19 in the two-dimensional 
quarter-filled band, to illustrate our point. 

For large U and V, the quarter-filled band is a CDW. The configuration 
with one extra particle is as shown in Fig. 9. This has configuration interaction 
with a higher energy configuration which is at energy 2V relative to it. If, 
however, one adds two particles which occupy neighboring empty sites in the 
CDW, the lowest energy configuration has configuration interaction with a 
configuration that is at an energy V relative to it. Thus, CI should give pairing 
in two dimension, though not in one dimension. This was actually numerically 
demonstrated by Mazumdar and Ramasesha 19 . The fallacy in this argument 
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becomes clear once three particles are added to the system as shown in Fig. 9 
Now, the lowest energy configuration has CI with a configuration at the same 
energy, so that a cluster of three particles is favored over a cluster of two and 
an isolated particle. Similar qualitative arguments favor phase segregation 
in related models. Thus, once again, thus, the pairing found within such 
calculations is a finite size effect. 

III. Conclusions 

The major conclusion of this presentation is that there are many dif- 
ferent finite size effects, and even when a given approach is satisfactory for 
the determination of one property of a given class of materials (e.g., the en- 
ergy of the 2A g in the infinite conjugated polyene), it does not necessarily 
mean that the same approach is satisfactory for other properties (in this case 
ground state bond alternation and the intensities of TPA). Finite size effects 
should be investigated for each different problem separately. Probably the 
most useful approach to avoid the pitfalls is to have a physical and intuitive 
mechanism of the effect in question in mind. However, this can hardly be 
prescribed as a general approach. One simple way to investigate finite size 
effects is to carefully investigate the size dependence of the property in ques- 
tion within the noninteracting model, for which calculations can be done for 
arbitrary sizes. Even though this procedure is simple (and an obvious step), 
finite correlated clusters are often compared to the infinite noninteracting sys- 
tem, simply because the latter is understood. Such a procedure can, and often 
does, erroneously ascribe effects due to finite size to correlations. 
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Table 



Table 1. The ratio of the negative and positive contributions to the 2A g TP A 
in long polyene chains within the extended Hubbard model. 
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Figure Caption 

Figure 1. The single-particle orbital electron occupancies of the lA g , \B Ul 
2A g and the 2B U . Only the highest (lowest) two valence (conduction) band 
levels are shown. All other valence (conduction) band levels are filled (empty). 

Figure 2. The intensity of the 2A g TPA in short finite chains within the Hiickel 
model. Here E g is the 1B U — lA g gap. Notice that the intensity increases with 
N. 

Figure 3. The intensity of the 2A g TPA as a function of N within the Hiickel 
model. TPA goes to zero at N — > oo. 

Figure 4. The ratio of the negative and positive contributions to the 2A g TPA 
as a function of N within the Hiickel model. 

Figure 5. The electronic energy difference per site between dimerized and the 
equal bond length even chains for the parameters of Eq. (6). The extropolated 
N — > oo values lies above the value for the dimerized infinite Hiickel ring. 

Figure 6. Pair binding energies (see Eq. (7)) for the (a) N = 4 and (b) N = 8 
rings. In each case the solid curve corresponds to the undistorted ring, the 
dashed curve to the Jahn- Teller distorted ring. 

Figure 7. The gain in electronic energy on Jahn- Teller distortion for the N = A 
ring, as a function of Hubbard U. The solid and dashed curves correspond to 
the 4 and 3 (5) electron cases respectively. 

Figure 8. Pair binding energies in the simple cube and three distorted forms. 

(i) , simple cube; (ii) - — , stretched cube, t\ = ti > £3; (iii) - 

— , Jahn- Teller distorted cube; (iv) — - — -, all three hopping integrals 
different. Notice absence of pairing at small U for cases (iii) and (iv). 
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Figure 9. Quarter-filled extended Hubbard CDW with (a) one extra particle, 
(b) two extra particles, (c) three extra particles (crosses are particles, dots 
empty sites). CT excitation in (a) costs energy 2V ', in (b) V, and zero in (c). 
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